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THE LAG-RANGIAN MULTIPLIER METHOD OF FINDING 
UPPER AND LOWER LIMITS TO CRITICAL STRESSES 
OF CLAMPED PLATES 

By Bernard Budiansky and Pal C. Hu 


SUMMARY 


The theory of Lagrangian multipliers is applied 
to the problem of finding both upper and lower limits to 
the true compressive buckling stress of a clamped rectan- 
gular plate. The upper and lower limits thus bracket 
the true stress, which cannot be exactly found by the 
differential-equation approach. The procedure for obtaining 
the upper limit, which is believed to be new, presents 
certain advantages over the classical Rayleigh-Ritz 
method of finding upper limits. The theory of the lower- 
limit procedure has been given by E. Trefftz but, in the 
present application, the method differs from that of 
Trefftz in a way that makes it inherently more quickly 
convergent. It Is expected that In other buckling 
problems and In some vibration problems the Lagrangian 
multiplier method of finding upper and lower limits may 
be advantageously applied to the ‘calculation of buckling 
stresses and natural frequencies. 


INTRODUCTION 


Many important problems that cannot be exactly 
solved by the differential-equation approach and must 
therefore be analyzed by approximate methods arise in 
the buckling and vibrations of thin plates. The theory 
of Lagrangian multipliers can be a powerful tool in the 
analysis of many of these problems. The present paper 
presents the details of application as well as the 
fundamental principles of the Lagrangian multiplier 
method by demonstrating the use of the method to obtain 
both upper and lower limits to the true compressive buck- 
ling stress of a rectangular plate clamped along all edges. 



2 


NAG A TN NO. 1103 


The procedure for obtaining the lower limit is 
similar to a method used by Trefftz (reference 1) and 
recently described by E. Reissner (reference 2). The 
present lower-limit method differs from that of- Trefftz, 
however, in a way that makes itr inherently more quickly 
convergent. The upper-limit procedure, which does not 
appear to have been presented previously, is simpler than 
the usual Rayleigh-Ritz method and may be expected to 
permit the computation of more accurate results with less 
labor . 

In a recent treatment of the problem of— compressive 
buckling of clamped plates, extensive calculations of 
lower limits were made by Levy (reference 3) hy means of 
a procedure equivalent to the Trefftz method. The result's 
were estimated by Levy to be within 0.1 percent of the true 
results. In order to illustrate the methods of the 
present paper, upper and lower limits to the buckling 
stress of a square plate are computed to within 0.1 per- 
cent of each other; a positive check on the accuracy of ' 
Levy's results is thus obtained. 


a 

b 

P 

t 

E 

D 

x 

7 

w 


SYMBOLS 


length of plate, in direction of stress 
width of plate, perpendicular to. stress 
aspect ratio (a/b) 

. thickness 
Poisson's ratio 
Young's modulus of elasticity 
plate stiffness in bending 


Et3 


12(l - tl 2) / 


plate coordinate in direction of stress 

plate coordinate, perpendicular to direc- 
tion of stress 

plate buckling deformation, normal to 
plane of the plate 
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lc 


V 

T 

a n 
b r 
a mn 

1, j, m, n, p, q. 
r, s 

5 mn 

\in 

(x } X y 


critical compressive stress, in x direc- 
tion 


critical compressive stress coefficient in 
the formula, cr^ = kf ■ ^ ■ j 

internal energy of deformation 
external work of applied stress 
Fourier coefficient of , cos nTr ^" 

Fourier coefficient of cos 


a 

rtry 

a 


Fourier coefficient of cos 222 cos 

a b 


even integers 
odd integers 

Kronecker delta (1 if m = n* 0 if m ^ n) 

_JL 

[(m 2 + n 2 {3 2 ) 2 - km 2 (3 2 ] (l + 6 m0 + 5 0n ) 

Lag rang! an multipliers 


THEORETICAL BACKGROUND 


Rayleigh-Ritz method .- The Rayleigh- Ritz energy 
method for determining the critical stress of a thin 
plate consists of the following steps: 

(1) The deflection surface of the "buckled plate is 
expressed in expanded form _as the sum of an Infinite 
set of functions having undetermined coefficients. In 
general, each term of the expansion must satisfy the 
geometrical boundary conditions of the problem. 

(2) The energy of the load-plate system is computed 
for this deflection surface and is then minimized with 
respect to the undetermined coefficients. 
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(3) This minimizing procedure leads to a set of 
linear homogeneous equations in the undetermined coef- 
ficients. These equations have nonvanishing solutions 
only if the determinant of their coefficients vanishes. 

The vanishing- of this stability determinant provides 
the equation that may be solved for the buckling stress. 

When the set of functions used is a complete set 
capable of representing the deflection, slope, and 
curvature of any possible plate deformation, the solu- 
tion obtained is, in principle, exact. Since, however, 
the exact stability determinant is usually infinite, a 
finite determinant yielding approximate results is 
used instead. 

Lagrangian multiplier method .- The Lagrangian 
multiplier method 'follow’s the general procedure outlined 
for the Rayleigh-Ritz method, with but one outstanding 
change. The restriction in step (1) that the boundary 
conditions be satisfied by every term of the expansion 
is discarded and is replaced by the condition that the 
expansion as a whole satisfy the boundary conditions. 

This condition is mathematically satisfied in step (2), 
during the minimization process, by the use of Lagrangian 
multipliers. 

The fundamental advantage of the Lagrangian multiplier 
method lies in the fact that, with the rejection of the 
necessity of the fulfillment of boundary conditions term 
by term, the choice of an expansion is much less restricted. 
In the clamped -pi at'e compression problem, a simple Fourier 
expansion may be used instead of the complicated func- 
tions assumed in the Rayleigh-Ritz analyses of this 
problem (references I 4 . and 5)» Furthermore, the orthogo- 
nality properties of the simple Fourier expansion lead 
to energy expressions of a simplicity that is Instrumental 
in permitting accurate computations . 

Approximate solutions of upper and lower limits . - 
The Lagrangian multiplier method, as well as the Rayleigh- 
Ritz method, gives a theoretically exact solution for 
the buckling stress; but ordinarily only approximate 
results are obtained because of the practical necessity 
of considering finite rather than infinite determinants . 

In the Rayleigh-Ritz method the approximate result Is 
always higher than the true buckling stress. In the 
Lagrangian multiplier method, however. It Is possible to 
obtain approximate solutions in two different ways. 



NAG A TN No. 1103 


5 


which permit the computation of a lower limit as well as 
an upper limit to the true buckling stress. As deter- 
minants of higher order are used to obtain approximations 
of higher order, both the upper-limit and lower-limit 
results approach the true buckling stress. Thus, the 
Lagrangian multiplier method can be used to provide a 
result to within any specified degree of accxxracy. It 
may be expected, furthermore, that a particular deter- 
minant in the Lagrangian multiplier method ought to 
yield a more accurate result than a determinant of equal 
order in the Rayleigh-Ritz method. 


LAGRANGIAN MULTIPLIERS 


The procedure used in applying the fundamental 
mathematical principles of Lagrangian multipliers is 
described in this section; a general proof of the 
validity of the method is given in the appendix. 

Let it be required to minimize a function of N 
variables - _ 

f ( x l> *3 • • • (1) 

where the x’s are not independent but are bound 
together by the relationship 



( 2 ) 


Lagrange’s method of simultaneously minimizing 
and satisfying the constraining relationship (2) is to 
minimize the function 


f - \cp 


with respect to the x’s. The quantity -K is the 
undetermined Lagrangian multiplier. The necessary con- 
ditions for minimizing f then become 
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£sf ^ dcp _ 
6x k 


( E- 1, 2, 3 . . N) 


cp = 0 


(equation (2 ) ) 


Note that these expressions are N + 1 equations in 
the N + 1 unknowns- x^, Xjj- and A. . 

If there are.. two relationships that constrain 
the x r s; that is, if 


°Pl( x l> x 2 * ‘ ‘ x It) “ 0 
q>,(xi, x 2 . . . x N ) = 0 


two Lagrangian multipliers are 
tion to be minimized becomes 

f - - 


then needed. 
X 2 cp 2 


The f unc - 


and the minimizing equations are 

dcp, dq?^ - - 

J2£_ - \ 1 — 1 - \ 2 . — £ =0 (K = 1, 2, 3 . . . N) 
6xj£ ^ x K 

<f> 1 - 0 
<P Q = 0 


The method is easily extended to cover the case of 
any number of constraining relationships. 


PRELIMINARY ILLUSTRATIVE EXAMPLE 

Before the main example Is given, a simpler buckling 
problem will be analyzed by the Lagrangian multiplier 
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method in order that the method of application of 
* Lagrangian multipliers may be most clearly presented 

without the obscuring details of analysis of more compli- 
cated problems. This elementary problem requires the 
use of but a single Lagrangian multiplier which leads 
to a single stability equation. _ 


7 



Consider a square plate, clamped along two opposite 
edges, simply supported along the other two edges, which 
is loaded in compression on the simply supported edges. 

(See sketch.) Prom the exact solution of this problem, 
the deflection surface of the plate is known to be sinus- 
oidal in the x direction. The deflection in the 
y direction, known to be symmetrical, must satisfy the 
clamped-redge boundary conditions; that is, zero deflection. 


w(x, 0) = w(x, a) = 0 

and zero slope 



0) 



a) = 0 


( 5 ) 

(k) 


The present method uses a cosine-series expansion, 
whereas the Trefft’z "procedure would use a sine-series 
expansion. The problem is solved by both methods for 
comparison. 
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Cosine -series solution .- In the cosine-series 
solution the expansion 

CO 

w = sin a n cos (5) 

n=0, 2,1|. ... a 


may represent the deflection surface having m half- •, 
waves in the x direction, since the Fourier series of 
even cosines is a complete symmetrical set. 

The boundary conditions (equation (i+)) on the slope 
are satisfied by each term of the expansion; however, in 
order that w satisfy the conditions of equation (3) on 
the edge deflection, it is necessary that 


00 



n— 0 ,2,1}. ... 


a 

n 


0 


( 6 ) 


Equation (6) is a constraining relationship on the a's 
and as such will be introduced in Lagrange's minimization 
process. 


As In the Rayleigh-Ritz method, the internal energy 
of deformation and the external work of the stresses are 
then calculated. Using the value for w as given by 
equation. ( 5 ) in the general formulas (reference 6, 
equations ( 199 ) and (201 ) (modified ) ) 


V = 


D 

2 




- 2(1 







X 



dy dx 


(7) 


T 


_ o*t r a r b 
“ —\ 


Jo 


,0 


(dZi 2 

\<W 


dy (he 


( 8 ) 
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gives 

■ v = ^f y (m2 + n2) 2 (! + 6 0n ) a n 2 
8a n=0,2,k ... 


mP-Tpo t \ / >2 

T = >_ (l + 5 0n ) a n ^ 

n— 0 , 2 , ... 

The usual Rayleigh-Ritz procedure requires that the 
expression 


V - T 


( 9 ) 


"be a minimum with respect to the a’ s . In the present 
example, however, the a' s are not independent hut are 
hound by equation (6). Hence, mathematically stated, 
the expression V - T must be a minimum subject to the 
constraint relationship on the .a’s 


CO 

/ an = ° 

n-0 , 2 , ip • * • 


(Equation (6)) 


Solving this minimization problem by Lagrange’s 
method makes it necessary to minimize 


co 



n — 0 , 2 , Ij. • • . 


( 10 ) 


with respect to the a’s. The necessary conditions Tor 
a minimum then become 
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/ w 

»' z 


6(V - T ) _ \n=0,2.1t . 


da 


3 


6a j 


= 0 (j = 0, 2, Ij. . . .) 


a n = 0 


n=0,2,lj. 


or, upon differentiation and simplification. 


(1 + 5 oj) [( m2 + J 2 ) 2 - 


m^k 


1+a 2 

aj ' K x= 0 


(11) 

( j = 0 , 2 , Lj. . . .) 1 


v~ 

z_ 

n=0,2,[j. . . 


a = 0 
n 


( 6 ) 


Solving equation (11) for aj "and substituting into 

equation (6) gives the stability equation that deter- 
mines k: ' 


X 

j-oT^rir 


'(m 2 + j 2 ) 2 - m 2 k] (l + 5 0j ) 


= 0 


( 12 ) 


Fob a particular number of half waves m, this 
equation may be solved by evaluating the series for 
several trial values of k and. interpolating to find 
the k that makes the series vanish. The correct value 
of m is that which gives the lowest value of k. For 
two half waves ' (m = 2) in the loaded direction, the 
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theoretically exact value of k = 7*69 (reference 6, 
page 345 ) I s obtained when only ten terms of equation (12) 
are computed. 

Sine-series solution (Trefftz method) .- The same 
problem will be treated in the manner suggested by 
E. Aeissner (reference 2), which is similar to Trefftz’ s 
method (reference 1). 


(13) 


The boundary conditions on deflection (equation (3)) 
are now satisfied term by term, but che conditions on the 
edge slopes (equation (ip)) are satisfied only by making 




w = sin 


mtoc 


b r ain ^ 


*=1,3,5 


f rb r = Q 

r=l,T75 ... 


( 34 ) 


Nova, the expression V - T is computed from 
formulas (7) and (3) by using the value of w given 
by equation (13); then-, by application of Lagrange's 
procedure , 


03 



r— 1, 3 ,5 * • • 


(15) 


must be minimized with respect to b s (s = 1, 3, 5, • * *) 
The Lagrangian multiplier is y. 


The minimization equations are 


Km 2 + s 2 ) 2 - m 2 k 


), q 2 

b D - ys = 0 ( s 


n^D 


1, 3, 5, 


r— 1 ,3,5 * • • 


rb r = 


. .)(16) 

(lip) 


o 
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Solving equation (l6) for b~ and substituting in 
equation (lip) gives as the stability equation 


CO 



3=1, 3, 5 




- m^k 


0 


(17) 


Comparison and discussion of results .- The series 
in equation ( 17 ) converges approximately' as l/s 2 , 
whereas the series in equation (12) converges approxi- 
mately as 1/jM-, Because of the more rapid convergence 
obtained in the stability equation, the Lagrangian 
multiplier method is preferably used to satisfy the zero- 
deflection condition rather than the zero-slope condition. 
Slope is the derivative of deflection, and/ in general, 
differentiation of a series makes it more slowly con- 
vergent. 


THE COMPRESSIVE BUCKLING OF A RECTANGUIAR PLATE 
CLAMPED ALONG ALL EDGES 


The previous elementary example required only a 
simple Fourier expansion and but one Lagrangian multiplier 
to satisfy the boundary conditions . The more difficult 
problem of finding the buckling stress of the rectangular 
plate clamped on all edges and loaded as shown In the 
accompanying sketch necessitates a double Fourier series, 
as weil as an infinite set of Lagrangian multipliers to 
satisfy the boundary conditions . 



* 
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Boundary conditions .- The boundary conditions of 
the problem are : 


Zero deflection, loaded edges 

w ( 0 , y ) = w(a,y) = 0 

Zero deflection, unloaded edges 
w(x, 0 ) = w(x,b ) = 0 
Zero slope, loaded edges 

~(0,y) = —(a,y) = 0 
ox ox 

Zero slope, unloaded edges 

dvr. bv, , N 

- — (x,0) = — (x,b) - 0 
cy Oy 


(18) 

(19) 

( 20 ) 


( 21 ) 


Four le r expans ions . - In order to achieve a rapidly 
convergent solution, the principles established by the 
preceding example are used, as the basis for ^choosing the 
Fourier expansions to satisfy, term by term, the condi- 
tions of zero slope .rather than "those of zero deflection. 

The buckling deformation corresponding to the lowest 
buckling stress is always symmetrical perpendicular to 
the direction of load but, depending on the aspect ratio 
of the plate, may be symmetrical or antlsymmetrical In 
the direction of load. Thus, for symmetrical buckling, 
let 


w= > y - (22)" 

m=0, 2,Lj- . . . n=0,2,ij. ... 
and, for antisymraetrical buckling, let 

CO - 00 

w= > "V . o M 2S«.HI (25) 

. a b 

r=l,2,5 «•• n= 0 , 2,4 ... 
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It is sufficient, for purposes of demonstration, 
to consider only the case of symmetrical buckling. 
Hereinafter, w therefore refers to the value given 
by equation (22), 

Energy expressions . - Using the expansion given by 
equation (22)' in the' evaluation of the general energy 
and work integrals of equations (7) and (8) gives 


Z 8aT' ^ S ( m2 + n2 ^ 2 ) 2 ( 1 + 6 m0 + 6 0n ) a mn 2 

m=0,2,i). ... n= 0 , 27 IT*»* 


T = 


rr^o^tb 

8a 




m ; 


— 0 ,2 , !(. , , , n— 0,2)1). , , , 


+ 6 On) 


OnJ %n 


Then 


v _ t = 2V 

8 a^ 


J 


m : 


I_ o 2 

nnn 


Am 


(2k) 


0,2, J4. » , , n— 0,2,1). *c. 


where 


1 


A, 


■mn 


(m 2 + n 2 p 2 ) 2 - kni 2 p 2 ] (l + 6 m0 + 6 Qn ) 


Note that V - T is independent of a nn , since 

— = 0. 

A oo 

Constraining relationships The boundary conditions 
of zero slope ( equations ( 20) 'and (21)) are satisfied 
by each term of the expansion of equation (22), but the 
conditions on deflection (equations (lo) and (19)) must 
be satisfied by the expansion as a whole. Substituting, w 
into equation ( l8 ) gives, along the loaded edges. 
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w(0,y) = w(a,y) 


J 


^0 + 003 . 

m=0,2,i{. . . 


V s T *-a 

m= 0 ,' 2 , 4 - • • • 


+ cos 


km 




m — 0 j 2 j 1_[_" ... 


a ml{. 


+ cos 


inz 


co 

I 




■+• • 


xn — 0 j 2. j Lk • • • 


= 0 


Since this Fourier series must vanish, each infinite 
series that constitutes a coefficient of a cosine term 
must vanish. (All. the Fourier coefficients of the 
Fourier expansion of the function zero are zero.) Hence 


U = 0, 2, I4. . . .) (25) 


CO 

^ = 0 
m= 0,27Ii. • . . 


By expressing the fact that there is zero deflection 
along the unloaded edges (equation (19)), it can be 
similarly shown that 




in 


= 0 


(i = 0, 2, I4. 


) (26)' 


n=0,2,!+ 


Now, V - T must be a minimum with respect to 
the a’s, which are bound by equations (25) and (26). 
As the problem now stands, however, it is not in the 
form to which Lagrange's minimization process can be 
applied since V - T does not contain &oo> whereas 
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the constraint relationships do contain age* Hence , 
ago is eliminated from the constraint relationships by- 

subtracting the first of -equations (26), the equation 
for i = 0 , from the first of equations (23), the equa- 
tion for j ~ 0 . The final set of necessary constraining 
relationships on the minimization of the energy expres- 
sion (2I4.) 


V - T 


Apb 

8 a 5 



> 


n= 0 , 2,4 


Aimi 


®mn 


then becomes 


^nO 

m=2,I|., 6 . . . 



n=2 j&S 

CO 

r 

m= 0,274 


\ 

/ 

n- 0 , 2,TJ.’ 


a 0n = 0 

axnj ~ 0 ( j = 2,k* 6 

a in = 0 ( 1 = 2, U, 6 


.)> ( 27 ) 


.) 


Theory of upper and lower limit solutions .- A 
theoretically 'exact solution to the problem would be 
obtained if the energy expression (equation (21^.)) were 
minimized with respect to all the a's and at the same 
time all the relationships (27) were satisfied. This 
result follows from the facts that: (a) the expansion 

of equation (22 ) is a complete symmetrical set, capable 
of representing the exact symmetrical buckling deforma- 
tion, and (b) the fulfillment of the conditions of 
equations (27) ensures that the boundary conditions are 
completely satisfied. An exact solution is not possible, 
however, because it would involve an infinite determinant, 
so that two different modifications of the ideal pro- 
cedure are used to obtain approximate results. One of 
these methods gives an upper limit .to the true buckling 
stress, whereas the other gives a lower limit. 
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An upper limit to the buckling stresa can be found 
by arbitrarily setting some a’s equal to zero, 
minimizing expression (2l|) with respect to the remaining 
a’s, and satisfying all the constraint relationships (27). 
An upper limit is obtained inasmuch as arbitarily setting 
some of the Fourier coefficients equal to zero has the 
effect of restraining the deflection of the plate, which 
in effect stiffens the interior of the plate and increases 
the stress required to buckle it. 

A lower limit to the buckling stress can be found by 
minimizing expression ( 2l \. ) with respect to all the a’s 
but satisfying only some of the constraining relation- 
ships (27). Neglecting some of the constraining rela- 
tionships has the effect of giving the plate greater 
freedom at the edges and hence reducing the stress 
required to buckle the plate . 

Lower limit solution .- In accordance with the require- 
ments for a lower limit,' the constraining relation- 
ships (27) will be satisfied only up to j = q and 
i = p. By Lagrange's minimization process, the function 
to be minimized is then 
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The a, ^2 • • * \i > M-g • * * i-h, are Lagrangian 

T. JH 

multipliers. The Equations for minimizing V - T with 
the constraining relationships ( 27 ) on the a ! s then 
become ■ . .. : 


- 22 - = 0 (in, n = 0, ' 2, i|. . . . ) 
6a mn 

Equations (27) 

Evaluation of CG/fta^ gives 


T 2 ! — %>n - “( 6 0n - 6 m0 ) ' ^ = 0 

1 ^ \m • 


(28) 


( 29 ) 


where appears only if 2 ^ n =; _q and ^ appears 

only if 2 "h. m p. From equation (29) 

7 a D+d , q+e = 0 ( d v e “ 2 > & • • •) 

A p+d,q+e * - 

Hence for any particular d, e, either 


or* 


A p+d,q+e 
ap+d,q+e = ® 


The first alternative, however, ordinarily would 
require k to be very high, corresponding to the 
buckling stress of a buckling mode with many waves in 
both directions. For the lowest buckling load, then. 


5+d,q+e 


- 0 


(d, e ±= 2, ip, 6 
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It is therefore necessary to be concerned v/ith only the 
other a’s, which, from equation ( 29 ), are 


a mn 


ha- 


^nO " 


rr^-Db 

- Ju? 


rAlft 


a 0 n 


_ _4a 


rAhb 


A mo( GC + hn) 

A 0 n(" a + X n) 


N 


> 


(n, n 7 ^ 0 ) 


(30) 


•In equations (30), does not appear if. n > q and p- m 

does not appear if m > p. 


Substituting the values of a's given by equation (30) 
back into the constraining relationships ( 27 ) up to 
j = q, i = p gives ,• 




a [ y V) + l + ^ pnO^m “ A 0n\ 

=2,u,6... n=2,I;,6.../ m=2,l+, 6... n=2,4,6... 


_A O0 a + / A mj^m + 




: §X6... 


A , = 0 
mj 

( j = 2 , 1;, 6. . 


A i0 a + Pi A in + ^3 A ir? v 'n ® 

n =0,2,1;... n=2,l;,6... 

(i = 2, k, 6 


(3D 


q) 


p) 


These equations form a set of 


homogeneous equations in a. 


P'2 • 


^(p + q) +i 

• • p'p } ^2 * * 


linear 




20 


NAG a TH No. 110$ 


Since when buckling occurs the a’s are not all zero, 
by equation (30), the Lagrangian multipliers .are not all 
zero. In order that equations (31) be compatible, the 
determinant of the coefficients of the Lagrangian 
multipliers must vanish* The vanishing of this stability 
determinant provides the de terminantal equation that may 
be solved for k by substitution of trial values and 
interpolation. 

That certain elements of the determinant consist 
of an infinite series of terms is evident; these 

series converge rapidly. Since such rapidly convergent 
series are calculable to any degree of accuracy, they 
may be considered as known quantities. Each value 
of A^ represents the potential-energy contribution 

of a term in the expansion for w; hence, the effects 
of infinite subsets of expansion terms enter into this 
solution. Thus, for p = q = 2, the expansion terms 
corresponding to the a's shown in figure 1 enter Into 
the solution; similarly, the terms represented in 
figure 2 enter into the solution, when p = q = Ij.. 

Upper-limit solution .- The lower-limit solution 
satisfied only some, of the constraining relationships ( 27 ) 
but assumed the existence of all the Fourier coefficients. 
If an upper limit is to be obtained, it will be necessary 
to satisfy all the constraining relationships while 
arbitrarily assuming some a's to be zero. 

As a direct result of the necessity of satisfying 
all the constraining relationships in the upper-limit 
solution, it’ "is found that the first of equations ( 27 ) Is 
redundant and may be discarded, since it is automatically 
satisfied when all the remaining equations ( 27 ) are 
satisfied. As a proof of this redundancy, the conditions 

CO 

*c — 

Aftij - 0 (j-2, Jp, 6.,.) 

CT1- 0727I+ • f • 


are summed over j and subtracted from the sum of the 
conditions 


cx> 



n=0,2,i).. . . 


( i = 2, ij., 6 . . , ) 
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over 1 and give 

CO m 


L in 




5 


0 £_ 




= 0 


i 2 ,]± , 6 < . . n-0, 2,lj. ... j —2 . . . m-0 ,2,1^ . . . 


Simplifying this equation 

QO CO 

'y' a i0 “ / a 0j ” 0 

}L\- > & « • • j -2 ,ij-, 6 » . . 

which is precisely the first of equations (27). 

It Is to he emphasized that the redundancy of a 
constraining relationship is a peculiarity of only the 
upper-limit solution, since, as shown by the proof given, 
the redundancy depends on the fact that all the con- 
straining relationships must be satisfied. 

With the elimination of the redundant condition, 
the necessary constraint relationships become 


sfcj = 0 (j = 2, It, 6 . . .) (32) 

m=0,2,4.. • ♦ 

CP 

a in " 0 (i = 2, k, 6 . . ♦) (35) 

n=0, 2,4*. . 


At this point, in accordance with upper-limit theory, 
it is necessary arbitrarily to set certain a ! s equal to 
zero. It Is possible to take advantage of the Lagrangian 
multiplier method by allowing infinite rather than finite 
sets of a’ s to exist and still to obtain a stability 
determinant of finite order. Thus, infinite strips of 
coefficients of the type shown in figures 1 and 2 can 
enter into the solution. In the lower-limit case, the 
existence of all coefficients was assumed, but the coef- 
ficients ap + £ q +e were proved. to be zero; in this 
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upper-limit solution, it will be arbitrarily assumed 
that these same a ! s are zero; thus 


a p+d,q+e 


= 0 


( d , e = 2 , k, 6 . . . ) 


The constraining relationships ( 32 ) and ( 33 ) become 

CO 

y a m j = 0 (j 3 2, k, 6 . . * q) 

m=0,2,I l , . « 

P 


/ ^3 

m=0,2,i+ . . « 


a in 


= 0 


n- « 

q 


a in = 0 


n=0,2,k. . . 


(1 = 2, ]+, 6 , . « p) 


(32a) 


= 0 . (j = q +-2, q + i|., . . .«). ( 32 b) 


( 53 a) 


(i = p + 2, p + 4, . . oo) (53b) 


The function to be minimized is 

CO CO 


q_ 


g = 


rr^Db \ 

8a 5 ' / ^ 

;u=0,2 ,4. . . n= 0 , 2 , k . . 


" a mn 


^ j a mj 

3 = 2 , k ^'6 • . . m=0, 2^q ... 


xa_ 


X 3 ^ &m 3 " ^ 


Pi 


l in 


j=q+2 , q+lj. . . . 111=0,2, i>. . . . 1=2, <4. , ^ . n=0,2,lf. , . . 


y H y a in 

i=p+2 , p+i; . . . n=0, 2 ,1 l . . . 
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The first double summation of this equation extends over 
onlv the values of m and n such that 


Setting 


6G- 


m = p if n > q 
n ^ q if m > p 


= 0 then gives, for all the a’s arbi^ 


0a mn 

trarily allowed to exist, 

*“» = + 

where \ Q and do not exist. Substituting back into 

the constraint equations (32a), ( 32 b), (33a), and (33b) gives 


\ 




ra=0,2,4.... m=2,4-, 6 . . . 


— 0 ( j — 2, 4-, 6 . • . q ) (5+a) 


A m< j + y = 0 U = 3 + 2, q +4- • . .w)(34.b^ 

m=0,2,4- . . . m=2,4, 6 . . . 




A in^n + ^i 


II— 2. j S • • • 

q 

r ™ ' 1 

^in^n + ^i 


A in = 0 (i = 2, 4-* 6 . . *p)(^a) 


n=0,2,4. . . 

__q 

'X A^ n = 0 (i = p + 2, p + 4 . ..co)(35b) 


n— 2,4, 6 • . « 


n=0,2,4-. . 
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These equations involve all the Lagrangian multi- 
pliers. They can be reduced -to a set of equations, 
however, in . . . A, q , t ^ ie i’oHowinS 

manner: - ^ 

Prom equation (35t>)j for i = p + 2, p + Ip ... 


A 4 X 
in n 


n=2 , 1. j. , 6 . ■ . 


X A in 
n=u,2,4.. . . 


Prom equation (3llb), for j = q + 2, q + Iq . 

P 

- / An j 

\ — ^ 9 9 9 


X M 

m=0 . . 


Substituting these expressions for and X_j into 

equations (3l}.a) and (35a), respectively, gives as the 
final stability equations: ^ ^ v 


^ ^ + y " y- A irij| 

171=0. 2, 1+. . . m=2, 1+, 6. . . m=p+2, p+i;, . . J 



a~2,"Lj, *6# • • 


<7277 


\ . 7, < \~ . f m=2,l+, 6. . < ! \ . 

/’ A in n " A in| p ~ j ^iy A in 

n=2, ii, 6 . . . n=q+ 2 >q+i 4 .. . . \ Amn J ri®0, 2,4. . . 

\in=d72, I*. . ./ 


p (j 2, 6 « . «q) 

\ 

X _ 

= 2jl+J 6. . < I . \ A _ A 


A in p + 7 A in 0 


(i = 2, k, 6 . . . p) 
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Equations ( 36 ) form a set of ij(p + q) linear 

homogeneous equations in ^2 . . . ^q, ^2 . . • ^p. 

The stability determinant is the determinant of the 
coefficients of the A’s and p’s. - 

It is of interest to note that in the usual 
Rayleigh-Ritz solutions only finite sets of expansion 
terms are ever taken into account, and tho order of the 
determinant obtained is ordinarily equal to the number of 
terms considered. It' is then reasonable that a particular 
determinant obtained by the Lagrangian multiplier method, 
which considers infinitely more expansion terms than a 
Rayleigh-Ritz determinant of equal order, may be expected 
to give a more accurate result. 


Numerical example .- For the case of a square plate, 
(3 = 1, upper and lower limits were computed. The 
results for the buckling-s tress coefficient k were: 


Approximation 

Lower limit 

Upper limit 

First; p = q = 2 

9-99 

10.11 

Second; p = q = !(. 

10.07 

10.08 


The expectation that the Lagrangian multiplier 
method should give closer upper limits than the Rayleigh- 
Ritz method, for a given-order determinant, can be con- 
firmed for this example. A second- order Lagrangia n 
multiplier determinant gives an upper limit of 
k = 10.11, whereas Maulbetsch (reference Ij.) and Smith 
(reference 5 ) use complicated deflection functions in 
the Rayleigh-Ritz method to derive third-order deter- 
minants that give, respectively, k = IO.I 4.8 and 
k = 10 . 11 . 


It Is seen that the second approximation, requiring 
the evaluation of a fourth-order determinant for the 
upper limit and a fifth-order determinant for the lower 
limit, definitely establishes the value of k to within 
0.1 percent. 

Levy (reference 3 ) used an Ingenious method of 
obtaining lower limits that is , In fact, equivalent to 
the Trefftz method of using double sine series and 
satisfying the zero edge-slope condition by the 
Lagrangian multiplier method. On the basis of com-_ 
putatlons involving determinants up to order twenty. 
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Levy concluded that his results obtained from tenth- 
order determinants are within 0.1 percent of the true 
results. Inasmuch as Levy obtained k = 10,07i|. for the 
square elate, the present relatively simple upper- and 
lower-limit calculations show that his estimated limit 
of error Is correct for this case. 


CONCLUSIONS 


1. The Lagrangian multiplier method can be used to 
compute accurate upper and lower.' limits to the compressive 
buckling stress of a clamped rectangular plate, thereby 
bracketing the true result. 

2. Prom a consideration of rapidity of convergence 
toward the exact solution in clamped-plate problems, It 
is preferable to use an expansion that satisfies the 
zero-slope boundary conditions term by term rather than 
the zero-def lection boundary conditions, 

3 . Because of the fact that the Lagrangian multiplier 
method permits the effects of infinite subsets of- expansion 
terms to enter into the solution, it is believed that a 
particular stability determinant derived by the Lagrangian 
multiplier method will, in general, yield a closer upper 
limit than that obtained from a determinant of equal 
order In the Rayleigh-RItz method, 

ip • It is expected that the method of Lagrangian 
multipliers may be useful- in the analysis of. .other sta- 
bility and vibration problems » In particular, the method 
may be immediately applied to the de termination of 
vibration frequencies of clamped plates, and to the 
determination of_buc kling stresses of clamped plates 
under compression In two directions. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 3* 
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APPENDIX 


GENERAL PROOF OF THE METHOD OF LAURA” GIAN MULTIPLIERS 


Let it be required to minimize 

f^x 1# x 2 , xj . . . xjf^ . ( A1 ) 


where the N x’s are. bound by the P independent 
relationships ( P < N) 


x 2, x 5 . . . X N ) ( J =' ■ 1 , 2, 3 ... P) 


(A2) 


It will be proved that the equations for determining 
the .minimizing values of the x’s are: 


££_ + *., + k 


OXi 


bx 


K 


2 6 


6<p p 6c?p 

+ . . • + X = 0 


■K 


P bx,. 


(A 3 ) 




(K = 1, 2, 3 . . . N) 


^j( x l» x 2' x 3 • • • x n) = 0 (elation (A2)) 

(J = 1 , 2 , 3 . . . P) 


The X* s are Lagrangian multipliers; these (N + P) 
equations determine N x’s and P Xi s . 

If the. values of only (N - P) x’s are known , 
the remaining P x’s are determined from the P rela- 
tionships (A2). For convenience, consider the last ■ 

P x f s in equation ( Al ) to be dependent upon the 
first (N - P) x’s. Then, for f to be a minimum" 
its first partial derivatives with respect to the 
independent x’s must vanish, or: 


\ 
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df df ^ x N-P-fl ^ x N-P+2 

+ + '■’ ■ 1 ■* + . . • 

dx M dxjv -p+1 dx M ^ X N~ p+2 ^x tr 

M 


Of 8 x n 
bx N bx M 


•(m = 1 , 2, 3 . , ; ( N - p)(Ai+) 


But each of these equations contains P quantities 
that cannot be directly evaluated - the derivatives of 
the dependent variables with respect to the independent 
variables. For each value of M , these P derivatives 
are determined by dif fersntiat.tng each of the P con- 
straint relationships (A2). with respect to x^, Thus, 

dc Pj 6x N-P+l 6cp j ^*N-P+2 

— ■ 11 “ ' ■ 1 1 " * * + —————— ■ L, ' r 1 ' 

dx M ~^ X N-P+2 6x M 


dcpj dx N 

+ . . • + =- 0 (J = 1 ; 2, 3 . . . P) (A5) 

dx N dx M 


Now, for each particular value of M, equation (AI 4 .) 
and the P equations (A5) make up a set of ( P ■+■ 1 ) 
linear homogeneous equations in the (P + 1) quantities 1, 


dx 


N-P+l 


‘'N-P+2 



Since these quantities 


dx M dx M bx M 

are surely not all zero, the determinant of— their coef- 
ficients nmst vanish. Hence * it is found that for f 
to be a minimum it must necessarily be true thatj 
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bf 

bf 

6f 

6f ! 



P+1 

a • 

4x N-P+2 

• , 1 

dx N ! 

j 



&c?l 


frq>i 


6x m 

6x N-P+1 

• a 

6x N-P+2 

9 

6x n 


dc? 2 

dc? 2 

6cp 2 

CbCp 2 

= 0 ( a6) 

dx M 

1 

^ X N-P+1 

a 

" 1 ""■""" a a 

^ x N-P+2 

a 

• " — 1 — 1 

(M = 1,2, 3 . . (N - P) ) 

a 

• 

dcpp 

• 

a 

6cp P 

< 

4 

bop 

a 

6cp P 


Ax m 

^ X N-P+1 

a a 

^ x N-P+2 

« "* ' ' ' 

dx N 



It will now be demonstrated that these necessary minimiza- 
tion equations will hold if equations (A3) hold. Inter- 
changing the rows and columns of the determinant in 
equation (Ao) gives: 


bf 

6cp x 

dcp 2 

6qp p . 


dx M 

6x M 

* X M ‘ * * 

6x m 


bf 

6c?i 

6cp 2 

6cp p 


^ X N-P+1 

6x N-P+1 

h < . — , . # # 

^ X N-P+1 

• 1 l_u 

^ X N-P+1 


bf 

6cp x 

dq> 2 

6q> P 

= 0 (A7) 

^ x N-P+2 

t 

^ x N-P+2 

4 

a • 

^ x N-P+2 

• 1 U " 

6x N-P+2 

(M =1,2,3 ...(N-P)) 

* 

bf 

* 

* 

dcpp 

• 

• 

6cp 2 



dx H 

dx H 

6xj, 

dx N 
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The vanishing of this determinant is, however. 


the condition of 

compatibility of 

the equations 

df 


+ + ... 

dcpp 

+ X p — - o 





df 

+ > 

+ \ 6cp 2 ♦ . . . 

- 0 

b x w _ P+1 

\ 

ox N-P+l 

° X N-Pfl 

r d x 

UX N-P+1 

df 

, 6< ?1 

. -v dc? 2 

dc()p 

M \ P .. B Pi 

^ x N-P+2 

° x N-P+2 

° X N-P+2 

/x# p w 

63C N-P+2 

* 

df 

dx K 

, 6c Pl 

+ /U — — 

bx N 

dcp 2 

+ Ap — + , . . 

Ox^ 

p 5T 

UX N 


precisely 




> (■*) 


when they are considered as linear homogeneous equations 
in the quantities 1, X^, X^, . . . X p . 

Since a determinant (A7) exists for each value 
of M up to (N - P), a set of equations (a 8) exists 
for each M. It is seen that In these sets only the 
first equation varies, since only the first equation 
depends upon M. Observation shows that— all the 
(N - P) determinants of equation (A7) can be derived 
from the set of N equations 


dcm 6cm 8cp 

— ±_ + X-j_ + Xg +.,'•+ = 0 (equation (A3)) 

d X K b*K bx K 6 *X 

(K a 1, 2, 3 ... H) 

by successively writing the determinants of compatibility 
of the last P equations with each of the first (N - P) 
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v equations in turn. It has thus been proven that if 

* equations (A3) are true, the minimizing equation (a6), 

eqtiivalent to equations (A'J), must hold. 

It is seen, however, that equations (A3) are 
N equations in (N + P) unknowns consisting of 
N x ' s and P X»s. The remaining necessary P equa- 
tions come from the original equations of constraint (A2) 
Hence, the simultaneous equations (A2) and(A3) 


df 

dx K 


, 6cpi 

+ A» ' ■ 

dx K 


+ X 


692 


+ X 


dc?p 


0 (equation (A3l ) 


(K = 1, 2, 3 . . . N) 



(equation (A2)) 

(j— 1, 2, 3 • • • 2) 


are necessary equations for the minimization of 
f rx]_, . . . x N \, which was to be proved. 
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Figure 1.- Four infinite strips of Fourier coefficients 

of expansion terms. - - ■ 
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